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Abstract: Using the AdS/CFT correspondence, we study the thermodynamic properties 
Q_i and the phase diagram of matter fields on (2+l)-dimensional defects coupled to a (3+1)- 

dimensional Af = 4 SYM "heat bath". Considering a background magnetic field, (net) 
quark density, defect "magnitude" 5N C and the mass of the matter, we study the defect 
psj contribution to the thermodynamic potentials and their first and second derivatives to map 

^ the phases and study their physical properties. 

We find some features that are qualitatively similar to other systems e.g. in (3+1) dimen- 
sions and a number of features that are particular to the defect nature, such as its magnetic 
properties, unexpected properties at T — > and finite density; and the finite 5N C effects, 
i e.g. a diverging susceptibility and vanishing density of states at small temperatures, a 

physically consistent negative heat capacity and new types of consistent phases. 
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1. Introduction 

In recent years, there has been a great amount of activity related to applying the AdS/CFT 
correspondence [1, 2, 3] to scenarios that may have some relevance in the context of ex- 
perimental physics. While the first and most common example, the duality between a 
stack of N c D3 branes generating an AdS§ x S 5 geometry in the decoupling limit and a 
thermal N = 4 SU(N C ) super- Yang-Mills theory on its boundary may be of limited ex- 
perimental relevance, for example fields transforming in the fundamental representation of 
the SU(N C ) have been studies by introducing a small number of Nf "probe" D7-branes 
into the background, covering all of the AdS directions [4, 5]. Certainly this setup is still 
significantly different from QCD, but this model and the T-dual D4-D6 setup have received 
great interest [6, 7], in particular the thermodynamics and phase structure [8] and the "me- 
son spectrum" [9] - hoping that some results obtained from AdS/CFT may be sufficiently 
generic, such that they also apply to QCD. 

More recently there has also been great effort on applying the AdS/CFT correspon- 
dence to condensed matter physics which may be more promising, since there is only one 
QCD but there are on the one hand many different strongly coupled effective field theories 
in condensed matter physics and on the other hand there exist a large number AdS string 
vacua [10]. The vast majority of these efforts have been related to 2+1 dimensional field 
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theories. The first example was an AdS^ x S geometry obtained from an M2-brane setup 
in [11] to study some transport properties. Since then, many interesting properties have 
been studied such as Hall conductivity [12, 13], superconductivity and superfluidity [14, 15] 
or the Nernst effect [16], mostly in setups that are not based on string theory and hence 
do not ensure that the systems are pathology-free. Another approach to move towards 
experimentally relevant field theories has been the construction of duals of non-relativistic 
CFTs [17]. A not yet satisfactorily addressed question is what the implications of the 
Fermi-Dirac distribution is in AdS/CFT [18]. 

The most common example of how systems with conformal symmetry arise in con- 
densed matter physics is the quantum critical phase. This phase arises in the context of 
a phase transition at zero temperature, a so-called "quantum critical phase transition" at 
a "quantum critical point". At this point, the system displays scale-invariant behavior as 
also in other phase transitions, but in contrast to phase transitions at finite temperature, 
it extends into a whole region in the phase diagram that may be described by a conformal 
field theory, the so-called quantum critical phase; shown in fig. 1. [19, 20]. 

All of the above-mentioned applications of the AdS / CFT correspondence to condensed 
matter physics have in common that they are described by 2+1 dimensional field theories. 
In our 3+1 dimensional world however, all 2+1 dimensional systems are strictly speaking 
defects. In some cases this fact may be less relevant and in other cases more relevant. Hence, 
it is interesting to study the physics of a 2+1 dimensional defect in order to explore what 
difference there is to purely 2+1 dimensional systems. Defect field theories are basically 
field theories in which matter that is confined to some hypersurface interacts via a field 
theory in the bulk. While there is some review literature in a soft condensed matter context 
(for a review see [21]) there seems to be not much review literature related to the defects 
and their aspects that we are interested in. Hence, a motivated guess may be that they 
have many properties in common with surfaces, which have been studied extensively. To 
illustrate their properties, we can look in fig. 1 at a generic surface phase diagram of some 
system described by a bulk coupling and a surface coupling J s - for a review on the 
subject see [22]. There we see that over most of the parameter range the surface and the 
bulk are in the same phase and display a simultaneous "ordinary" phase transition. As 
we tune the surface coupling beyond a "special point", which is some critical multiple of 
the bulk coupling, the phase transitions on the surface and in the bulk separate into a 
surface phase transition and an "extraordinary" phase transition in which there is a phase 
transition only in the bulk. It is obvious from the ratio J s /Jb m this regime that the 
ordered phase on the surface extends to higher temperatures than the ordered phase in 
the bulk. However it is quite interesting that this splitting of phase transitions typically 
occurs as J s becomes greater than Jb and hence there is no "mirror symmetric" version of 
this plot. 

In this paper, we study the thermodynamics and phase diagram of the defect. This 
means, we have to consider the contributions to the thermodynamic and statistical quan- 
tities of the full (defect +bulk) field theory that are extrinsic in the area of the defect. 
Certainly, we will miss contributions that are extrinsic in the volume of the SYM, but 
these depend only on a "topological" parameter that is related to a difference in the level 
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Figure 1: Left: A generic quantum critical phase transition and quantum critical phase. Right: 
A generic surface phase diagram. 



of the gauge group N c — > 5N C between both sides of the defect. As indicated in [23], 
this parameter is also related to a "width" of the region around the defect. There it was 
observed from the quasiparticle spectrum, and the length scales found in various regimes, 
whereas here we will discuss it from a more direct point of view. It will also be relevant 
when we argue about the positivity of the entropy and the heat capacity. As in [24] , the 
other parameters are the net baryon number density p, an externally applied magnetic field 
B and the "quark" mass M q . 

To study the thermodynamics, we obtain the free energy from the euclidean action, 
from which we obtain the chemical potentials and response functions. We also demonstrate 
how to obtain the second derivatives of the thermodynamic potential in closed form, i.e. 
without having to use numerical derivatives; even in the case of finite M q . 

Using those quantities, we will map out the phase diagram and point out some impor- 
tant differences arising from the choice of the "free" thermodynamic variables. We then 
argue that the blackhole phase studied in [24] is indeed the thermodynamically preferred 
one and study the physical properties of the different phases. We will then also find in some 
region of phase space a first-order transition between two distinct regions of the blackhole 
phase. It will turn out that this is in some sense a smooth continuation of the blackhole- 
Minkowski embedding phase transition, that was extensively studied in 3+1 dimensional 
systems and beyond [8, 7, 25, 26]. In contrast to the observations in those systems, it will 
turn out that a finite density Minkowski embedding can be realized in our case also at finite 
density, even though it is only metastable. Above the phase transition, the physics that we 
find will be dominated by some simple scaling laws, but below, in the small-temperature 
regime, we will find some surprising non-trivial effects. 

The gravitational setup of the defect is a stack of Nf probe D5- or D7-branes inserted 
into the background of a stack of N c D3-branes and is well known from the literature 
[4, 27, 28, 29, 30]. As in [23, 24], the difference SN C in the level of the gauge group of 
the 3+1 SYM will be introduced by an additional flux on the probe brane in the compact 
sphere, which also stabilizes the D7 setup. Similarly, the finite magnetic field and net 
density are introduced using the well-known duals of a magnetic field and an electric field, 
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respectively, in the world volume of the probe brane. The finite quark mass will be obtained 
by a deformation of the embedding in the compact sphere in the same fashion in which it 
was done in the duals for 3+1 dimensional QCD-like systems [8, 25, 31, 7]. 

The (3+l)-dimensional system of the N = 2 gauge theory constructed with parallel 
D7- and D3-branes (see e.g. [7, 6, 8]) is a generically good reference to which to compare 
our analysis. There, it was found that if a finite M q is introduced, the scale M iun ~ M q j\f\ 
plays a special role in this theory. At vanishing density, the picture is as follows: First, the 
"mesons" , bound states of a fundamental and an anti- fundamental field, are deeply bound 
with their spectrum of masses characterized by M fun [9]. Next at a temperature T ~ M fun , 
the system undergoes a phase transition characterized by the dissociation of the mesonic 
bound states [8, 7]. In the presence of background magnetic fields, it was found that the 
phase transition moves to larger Temperatures [32]. Some aspects of the (2+1) dimensional 
conformal matter in the presence of a magnetic field have also been studied in [33], which 
appeared on the arXiv during the closing stages of this work. 

The outline of this paper is as follows: In section 2, we briefly review the construction 
of the defect setup and the AdS/CFT dictionary for the background quantities and point 
out several problems that arise in the D3-D7 setup - which motivate us not to pursue 
the massive D7 case. We then discuss the thermodynamics in section 3, starting with a 
discussion of the thermodynamic potential and the thermodynamic variables in section 3.1, 
from which we then obtain the response functions in section 3.2, which we discuss first in 
the massless case in 3.3 and then in the massive case in 3.4. In the latter case, we first 
outline the phase diagram in section 3.4.1 and then discuss the physics of the response 
functions in 3.4.2. 

2. Gravitational setup 

In this section, we briefly review the string theory setup of the defect that was outlined 
in detail in [24]. We also point out some minor differences in the case of non-blackhole 
embeddings. 

2.1 D3 = (N=4 SYM) background 

The well-known AdS§ x S 5 background from the decoupling limit of N c D3 branes at finite 
temperature T = can be written as 

ds 2 = J (-Hr)df + d4) + £ + r 2 ^) , Qg, = (2.1) 

and corresponds to an N = 4 SYM theory on the boundary with U (N c ) gauge group [1, 2, 3]. 
We work, as usual, in the limit N c — > oo, Yang-Mills coupling g\ M = 2irg s — > and t'Hoft 
coupling A = gy M N c — > oo, i.e. we are in the supergravity limit L 4 = Ang s N c l 4 — > oo. 
Considering only T > allows us to go to dimensionless coordinates u = ^, t = j4, x = 
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The fields in the SYM all transform in the adjoint representation of the SU(N C ), 
and we follow the well-known probe brane approach of introducing Nf families of fields 
transforming in the fundamental representation, in the "quenched approximation" Nf <C 
7V C . Commonly, one inserts for example Nf D7 branes in the D3 background parallel to the 
D3's [5]. Generalizations to defect configurations have been considered in the literature 
[4, 27, 28]. In the latter cases, the intersection overs only part of the flat directions, 
creating the defect field theory, where the fundamental fields are only supported on a 
subspace within the four-dimensional spacetime of the gauge theory. In the case of our 
(2 + l)-dimensional defect, this can be done with the following configuration, 
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The D5-brane construction reduces the supersymmetry from N = 4 of the D3 back- 
ground to J\f = 2 and the dual field theory is now the SYM gauge theory coupled to N t 
fundamental hypermultiplets, which are confined to a (2+l)-dimensional defect. In the D7 
case, all supersymmetry is broken and the defect CFT contains only N f flavors of fermions, 
again in the fundamental representation [29] . It turns out that the lack of supersymmetry 
in the latter case will manifest itself with the appearance of instabilities. In [23], we showed 
how this can be avoided, and we discussed in [24] how this instability becomes apparent 
in the scaling dimension of the scalar field that corresponds to the deformation of the S A 
of the D7-worldvolume inside the 5 5 background. There, we also discussed some problems 
related to the reliability of the quenched approximation that we consider in this paper. 

In this limit, N t <C N c , the D5-branes may be treated as probes in the supergravity 
background, i.e., we may ignore the gravitational back-reaction of the branes. For the 
D7-branes, however, this is only true locally and not in the asymptotic regime. 

2.2 Introducing the defect 

Considering only the U(l) subgroup of the U(Nf), the action of the D5 brane in a D3 
background is just the DBI action plus a Chern-Simons term 

S = -T 5 N f [ y/-det{P[G] + 2irP s F) + T 5 N f [ C (4) A 2irl*F , (2.4) 

J D5 JD5 

where the factors of Nf arise from taking the trace over the flavor degrees of freedom. 

Assuming translational (in the flat directions) and rotational (on the sphere) symmetry, 
the induced metric of the embedding (2.3) can be written in the form 



ds> = * (-(i - u* )d p + dii + (i + (i - </) (zx'iuf + - 2 r^)) 

+ u 2 {l-^>{u) 2 )dQl) , (2.5) 
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where we used the inverse radius u = ^ and \ = t~? = ttTz. Here and in the rest of the 
paper, we also use the notation ~ to denote quantities that are made dimensionless with 
appropriate factors of ttT unless explicitly noted otherwise. 

Now, we wish to consider a few interesting background quantities. Firstly, we can turn 
on a magnetic flux on the sphere, F = -^df^, which corresponds to inducing an extra 

number of colors, SN C = q G Z, parametrized by the quasi-continuous / := jjfj-q = , 
on one side of the defect and causes the probe brane to bend towards that side. Both from 
the embedding geometry d u x{u) = , ~f , and from the resulting quasiparticle spectrum 

in the field theory, it was argued in [23] that this flux also introduced a finite width, Az, 
of the defect. 

Using the AdS / CFT dictionary in [3] in analogy with the 3 + 1 dimensional system, 
e.g. [25, 32, 12], we find the gravity dual of the baryon density 

5S V\N C T 2 ~ 

\mid u Ao{u) (2.6) 



5A b Q dy ' 47T u->0 

and magnetic field B to be related to a non-trivial U(l) background on the brane: 

F\ u ^ = -E{u)dt Adr + Bdx A dy =: F E + F B . (2.7) 

We can also define the (asymmetric) background metric 

G = g + F. (2.8) 

A finite quark mass M q and dual condensate C can be associated with a non-trivial 
embedding *(n) by analogy with the 3 + 1 dimensional D3 — D7 system [8, 25, 31, 7]. This 
condensate has on the one hand an interpretation as a chemical potential for M q and on 
the other hand is considered in QCD contexts as the order parameter of chiral symmetry 
breaking. Parametrizing the S 5 as d£l 2 = dip 2 + cos 2 tp dQ, 2 , + sin 2 tp dQ, 2 . and putting the 
D5 on the first S 2 , such that simp =: the DBI-CS action becomes 

S = 4irL 2 T 5 J da 4 (y/-detGy/(l - * 2 ) 2 + f 2 + fu A d u z^j (2.9) 

= 4ttL 2 T 5 j da 4 (V^detG^/l + F 2 ^Jl + F| v / (1 - * 2 ) 2 + f 2 + fu^z^ . 

The background solution is then 



d u A (u 



B = const. (2.10) 

pyi+TV 1 - ^ 2 ( u ) + u 2 h(u)V{u) 2 



y/l - t>{u) 2 ^l + (/ 2 + (p 2 + B 2 ){1 + / 2 ))u 4 + (1 + B 2 u 4 )^{u) 2 (^{u) 2 - 2) 



-fJl - ^ 2 (u) + u 2 h(u)W(u) 2 

OuX ~ 



vrT^l - ^{u) 2 Jl + (f 2 + (p 2 + B 2 ){1 + f 2 ))u A + (1 + B 2 u A )^{u) 2 {^{u) 2 - 2) 
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where h(u) = 1 — u 4 , p := 7= 4?rp „ n and -B = ? ^L vl . Now, the width of the defect from the 
brane picture, z ma:r := lim u _>i z(u), decreases with increasing density or magnetic field, 
in agreement with the quasiparticle spectrum observed in [24]. 
The equation of motion for ^(u) becomes 

2(1 + BV)(1 - ^ 2 ) 3 + u 2 {\ - u 4 )(l - {f + {p 2 + B 2 )(l + f 2 ))u 4 + (1 + B 2 u 4 )^ 2 {^ 2 - 2))V 2 



u 1 



(l-^ 2 ) J(l-* 2 )(1-^ 2 + {v?- u e )V 2 ) (1 + (/ 2 + {p 2 + B 2 ){1 + / 2 ))u 4 + (1 + B 2 u 4 )^ 2 (^ 2 - 2)) 



■ ,, T , 1-^ 4 / l + (/ 2 + (p 2 + ^ 2 )(l + / 2 )K + (l + W)^ 2 (^-2) 

"" 1 u 2 V (i-* 2 )(i-^ 2 + (u 2 -u 6 )^' 2 ) 1 

which has no analytical solution, except for some limiting cases. For u — > 0, it is easy to 
see that the solution becomes 

* ~ fhu + cu 2 , (2.12) 

where rh and c are dimensionless free parameters that are determined by the boundary 
conditions at u = 1. Following arguments of the T-dual (3 + 1) dimensional D3-D7 [8, 25, 
31, 7], the quark mass M q and condensate C are given by 

M i = = ^2^^ ^ ° = ^^ r l N A T ^ = ^T 2 NfN c c . (2.13) 

In section 3.1, we discuss this more in detail and verify that C is indeed the dual chemical 
potential to the mass. 

In order to find the solution for the full geometry for a given mass however, we need 
consider the equation near the horizon, where (2.11) reduces to first order, 

VI = l - (l + g 2 )*o(l-*g) 2 (214) 

2(l_* 2 )2 + / 2 + p2 + £2( 1 + (1 _^2 )2 ) ' V- ) 

effectively relating rh and c. Hence, the only remaining boundary condition at the horizon 
is ^m^i = ^o- We then have to find recursively V&o for a given value of rh. 

So far, this embedding is the same as the the one discussed in [24]. Now, however, we 
also want to consider non-blackhole embeddings. In this case, the probe brane does not 
extend to the horizon u = 1, but closes off, i.e. VP — > 1, at some finite value of u, u max , 
above the horizon. The above discussion carries over, but now we are tuning u max instead 
of ^o- Just as at the horizon, the equations of motion reduce to first order at this point, 
and the boundary conditions at u rnax can be read off from the expansion 

lT , _ -, Umax — U {Umax — u) 2 V? max (3 + 13lft max ) 3 

~ U Jl+^ ) + 6(1 -ifi + I 2 + U ^ U max-U) . (Z.Lb) 

While in the case of the black hole embedding, the boundary condition at u = 1 and the 
equations of motion ensure that VP £ [0, 1[ over the whole range of u, we now have to discard 
unphysical embeddings in which VP < at some values of u. These solutions appear at 
large magnetic fields and are in practice related to having ^'\ u _ ¥0 < 0, i.e. they would 
correspond to a negative value of the mass parameter. 
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In the non-supersymmetric case, the magnetic flux on the compact S 2 is replaced by 
an instanton on the compact S 4 , and the coupling to the five- form flux comes now via the 
Chern-Simons term 

(2vr^ 2 



2 



T 7 N f / C (4) AF AF . (2.16) 



D7 



We also parametrize the S 5 slightly differently as dfl 2 , = dip 2 + cos 2 ip dQ 2 , such that the 
induced metric is 



ds 2 = §(-(!- t^d? + dxl + (l + (1 - u 4 ) (,» 2 + x _ ^ (u)2y ; l _ u4 

+ u 2 (l - ^(u) 2 )dnl) . (2.17) 



-T 7 N f [ F^AAAF -»• 8T 7 iV f ^ / du ; £ / y4 A F . (2.18) 



It turns out that a qualitative difference arises only in the massive case, as in that case the 
five-form pulls back to the S 4 factor of the brane as 

2 

This was used in [30] to obtain a Hall effect in a setup that would be considered from the 

perspective of this paper unstable. 

It is a straightforward exercise to verify that it turns out that in the massless case, the 

embedding and the solutions for background fields take precisely the same form as the ones 

in the supersymmetric case, under an appropriate re-definition of the flux parameter f\ = 

, 2 in terms of the instanton number = q 7 = ^ <L 4 TrF AF G Z. Furthermore, 

y/l+2\Q\ 67r y 8?r Js ' 

even the value of the action, which reads now <f s4 d 4 £l = (A r f L 4 (l — ^ 2 ) + 6n 2 £ A \ q 7 \ ) , is 
identical, modulo an overall factor. 

A problem related to the stability of the D3-D7 solutions was pointed out in [24]. The 
"smooth" instanton solution discussed in [34, 35] that preserves the rotational symmetries 

of the S" 4 , and hence also the symmetries of our field theory, limits — q — - > q 7 . As 
discussed in [23], the scalar mode of the supersymmetric setup is always stable, but the 
non-supersymmetric setup requires f 2 > 49/32 for its mass to be above the BF [36] bound. 
This stabilization becomes inconsistent, however, in the light of backreaction, as the D7 
brane cause an asymptotic deficit angle of Nf/12, and we would require Nf 3> 1 for a 
finite value of fa. Hence, we do not follow the path of the massive D3-D7 embeddings 
that are distinct from the D3-D5 case. Just for curiosity, this instability is reflected in the 
asymptotic solution for \&, which becomes in the non-supersymmetric case \l/(u) ~ u a± 

with a± = | ± y/ 4Q 2+4( 7^ ~ raising interesting questions about the nature of the dual 
"mass" operator. Here, we see that satisfying the BF bound for \E' precisely corresponds 
to real values of a±, i.e. to a non-oscillatory solution, and the unstable solution would 
be non-physical. In [30], it is, however, given some interpretation in the context of the 
quantum Hall effect. 

3. Thermodynamics of the defect 

Now, we are ready to compute the contribution of the defect to the extrinsic thermodynamic 
quantities that are localized on and around the defect, defined uniquely such that we assume 
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the bulk contribution to be translationally invariant anywhere else than at the defect, where 
it is at most discontinuous, and the defect contribution to vanish far away from the defect. 
This allows us to study the phase diagram of the matter, the level of stability of the different 
phases, and their physical properties. 

Naturally, we will implicitly consider the density of these quantities per unit area of 
the defect in terms of the boundary metric, which removes the divergence from the infinite 
volume factor of the integral. We will miss however any contribution to possible changes 
to the asymptotic characteristics of the 3+1 SYM. Since the asymptotics however should 
only depend on the topological properties of the defect, i.e. the flux parameter / and not 
the local properties of the embedding, this is still sufficient to map out the phase diagram 
and discuss its properties. 

As pointed out in section 2.2, our discussion of the massive embeddings is limited to 
the D5 case. 

3.1 Free Energy and Thermodynamic Variables 

As a starting point, we can straightforwardly compute the free energy via the standard 
procedure from the Euclidean action, I e , [8] using 

F = TI e , I e = f ^ £ e + hdy , (3.1) 

where the boundary terms 

hdy. = -\^f+l^ 2 (3-2) 

are dictated to us by consistency [37]. The other boundary terms of [37] do not contribute 
in our case because of isotropy in the flat directions and rotational symmetry on the 
sphere. Since the determinant of the boundary metric, 7, vanishes on the horizon, only the 
asymptotic boundary contributes. In the case of non-blackhole (Minkowski) embeddings 
there is only one (i.e. the asymptotic) boundary, provided the action is consistent at u m i n 
and the point at u m i n is included, as will be discussed below. 

In the next step, one can construct other thermodynamic quantities, such as the en- 
tropy, S = — the energy E = —F + TS or the heat capacity cy = §f;|y Since the 
defect is in thermal equilibrium with the bulk, and also the extrinsic curvature of the hori- 
zon in the brane geometry is the same as in the bulk theory, the choice for the temperature 
is obviously the bulk temperature T = ^gj. There is a slight ambiguity as to what one 
considers to be the thermodynamic volume V . One could either consider the defect as an 
isolated thermodynamic system, embedded in the SYM heat bath, or as part of an overall 
system. In the former case the volume is either the 2-dimensional volume J dx dy with the 
effective width of the defect considered to be an "internal" degree of freedom or alterna- 
tively the 3-dimensional volume J dxdydz over a finite width Az, e.g. Az = J z'(u)du, 
with z'(u) given in (2.10). In the latter case, however, one considers a large volume of 3d 
SYM plus the defect, and the quantities that we are studying are just the contributions 
that are extrinsic in the two dimensions of the defect, and independent of the extension of 
the volume in the z-direction - in the limit of placing the boundary of the volume far away 
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from the defect. To make this more explicit, we can look at the variation of the euclidean 
action, which gives 

or in terms of field theory quantities 

SF = 2C5M q + P 5fi + F f 5z max . (3.4) 

In principle we might expect a term oc mSm that would be divergent. This term however 
cancels because of the renormalization. Also the variation of the magnetic field does not 
contribute because of its tensor structure. Ff is just defined as the change of the 3+1 free 
energy density F$ym due to varying N c — > N c + q, in dimensionful units: 

Ff = = -<^N c N f T* = — ^ fN c Nf VXT 4 . (3.5) 

If / ~ 0(1), then term is suppressed by a factor A -1 / 2 compared to F$ym- Noting the 
fact that the SYM background is isotropic such that the pressure equals the free energy, 
Fsym = PsYM-, we could interpret FfS Az as a work term for the case of the isolated 
defect. This demonstrates nicely that this case is inconsistent, since want to consider 
the defect system "on-shell" and study thermodynamic processes obviously with a fixed 
gauge group, i.e. at constant /. This implies that we cannot use z max as an independent 
thermodynamic variable. Hence, we need to do a change of variables in the thermodynamic 
potentials, corresponding to a Legendre transformation of the action. Since we also want 
to consider processes at fixed baryon density on the defect, rather than at fixed chemical 
potential, as in [25], we do Legendre transformations in \ an d A t : 

C e = C e + pA' t - fx' or I e = I e + pA t - fx . (3.6) 

Now, the variation of the free energy is 

SF = 2C5M q + p5p + z max SF f , (3.7) 

which implies F = F(M q , p, f,T). Obviously we consider / fixed, even though it is not 
inconceivable to have processes in condensed matter physics that change the effective gauge 
group. 

In the case of Minkowski embeddings, the significance of the Legendre transformation 
can be seen nicely from the brane tension in the z direction at the endpoint of the brane 
Umax that can be straightforwardly computed from 

Tz = 9zz—, = QzzZ — x , (3.8) 

og Z z og rr 

where g is here the D3 background metric. This can be shown upon substitution of (2.10) 

and the boundary condition ^ ~ 1 — (u max — u) ,-}. 4 — ^ to match precisely the tension 

of qNf D3 branes. 
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In particular, we can look at the source term that corresponds to attaching an appro- 
priate stack of D3 branes in the flat directions at the endpoint of the probe branes and 
balances this tension to allow for a consistent static embedding, 

C sourc lu) = -f z{u)5(u max - u) . (3.9) 

After integration in the radial direction, this term gives precisely the same contribution to 
the action as the term — fx' that we added for the Legendre transform. An interesting 
comment to add is that in this Minkowski embedding, the radial tension r r vanishes at 
the endpoint, and is entirely generated by the 5-form flux acting on the brane. Thus, 
it is actually possible in this defect setup to construct an at least metastable Minkowski 
embedding at finite baryon density within string theory - in contrast to the setups of 
fundamental matter in (3+1) dimensions in [25]. 

3.2 Dual potentials and response functions 

As some physical quantities of relevance, we will obtain the entropy S = — fy|y, the 



total energy E = F + TS, heat capacity cy = §f?|y = — §jf 



potential fi = ^ and magnetization M = §g| T y 

The latter two quantities are straightforward, since we only need to keep in mind the 
temperature scaling and normalization, such that 

dF 4vr dF , , 

M = g£ = ' ?Z = : ^ N f' T M. (3.11) 
dB (7rT) 2 dB tt 2 V ' 

In principle, there would also be a contribution from the variation of the embedding ty{u) 
via a term &d Pt B^, however, it turns out that it does not contribute for the following 
reason: Since we did a Legendre transform in (A' t ,p) and (x'>/)> the variation ^ at 
constant /, B and p is on-shell, i.e. only a boundary term contributes: 



, baryon number chemical 



dF I 



^1 + ?h*L \ 



(3.12) 



My. 



In principle, this term depends on rh and c. Keeping rh fixed, c will generically depend on 
B and p and hence 5^ = (^d~gcj u 2 5(p,B). However, it turns out that in an expansion 

around u = and ignoring overall factors, we have + 9 q^, v ~ ^ +m 2 , such that any such 
term will not contribute in the limit u — > 0. For black hole embeddings, there is obviously 
also the second boundary at the horizon, but this contribution vanishes since the boundary 
metric vanishes, 7 — >■ 0, as u — > 1. In the case of Minkowski embeddings, there exists only 
one boundary (the asymptotic one), as we learned above that the Legendre-transformed 
action is consistent at the endpoint u max and corresponds to including this point in the 
integral. To demonstrate that there is indeed no contribution from the region at Urnax ? wo 
can use the expansion (2.15) of \I> near u max . This relates any change u max — > u max + 8u r , 



'"max 
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to some 5^. If we then were to exclude the point at u max and evaluate the integral only 
up to some u t = u max — e, the contribution from any on-shell variation of the scalar ty, 
— f§7 \ u=u S^(u e ) cancels the corresponding extra contribution to the integral, £ e (it e ) ^t^j 
as we take e — > 0. 

Computing the entropy requires a few more steps. In order to compute the temperature 
derivative, we proceed by computing the integral (3.1) in terms of the dimensionless versions 
of the coordinates and fields. We then consider only the indirect temperature dependence 
of the terms in the integral and an overall explicit temperature factor which we find by 
dimensional analysis to be T 3 . Keeping p, B, M q and / fixed, the temperature dependencies 



that we will need are given by ^ 



9 p dB 

6 T , Qrp 



2f and f|| M = -f . The variation 



B 

of the action with respect to p and B are straightforward and defined in (3.10) and (3.11), 
and the variation with respect to ^ gives on-shell 

S y6i e ac e „. . r„ „. . 0i Mu „. 2 2 / 3 



5V 



■» + /*» + -£» - -h^ m - <313) 



where E-q, is the equation of motion for Hence, in the limit u m i n — > 0, only the linear 
term in 5^f contributes and we do not have to worry about the non-trivial temperature 
dependence of the condensate. Finally, we have to worry about the temperature dependence 
of the boundary of the integral. The boundary at the horizon is fixed in terms of the 
dimensionless coordinate at u = 1, so there is no contribution from the horizon. For 
the asymptotic boundary, we can compute, in dimensionless variables and ignoring overall 
factors: 



dl. 



du r , 



u=u„ 



£ e + d Umin I bdy . = -^ + ^^ + B 2 + f + c 2 -^+0(n mm ) (3.14) 



Since we take the limit u m i n —> 0, any temperature dependence at fixed maximum dimen- 
sionful radius r max is proportional to dru m i n = Ur J^ n and hence, there is no contribution 
from this boundary either. Putting all the non-vanishing contributions together, we arrive 
with 

3F 2 fdF „ dF ~\ 5I £t .„ 1 . 

S = —=r - ^ t^pH + -v^^bdv or (3.15) 

T T\dp dB J ^ v ' 

TS = -3F - 2(ip - 2MB + M q C . (3.16) 

Finally, we can compute the total energy, E = F + TS and the heat capacity 

dE 



CV = df 



rr, d 2 F 
T 



v 



(3.17) 



v 



The second temperature derivative is computed in the same way as the first derivative for 
the entropy. The most straightforward way is to take the temperature derivative of the 
entropy. Keeping in mind that the total temperature derivative of the implicitly present 
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equations of motion vanishes, the remaining terms are: 



d 2 F 
dT 2= 



2/3. 



2B 



&r + fd- p + —d § + 



d^ 5 
~d~T6^ 



... - 2p„ 2B 
T + y 9 p + ~Y d B 



d M q C 



(3.18) 



F — 3- 



M n C 



J>2 



+ 



— ^ + — ^ 




where the derivatives are considered to act on anything towards their right. Here, we find 
that we cannot avoid computing the temperature dependence of the embedding \£ and the 
condensate. This can be done in closed form, without resorting to numerical derivatives. 
To do so, we expand the temperature derivative of the equations of motion for ^ to first 
order in temperature, obviously keeping the appropriate dimensionful quantities fixed. The 
term linear in temperature then gives us an inhomogeneous linear second order equation 



for ^ 
iUi dT 



~ 5 E-q, 5 Ety 
B — — + p- 



8B 



5p 



T ( 5 E<$, d~§> 5E^d^' 5Eyd^" 
JWdT + JW~dT + JWdT 



(3.19) 



and a boundary on the horizon condition fixing (j^) 1 d u (^fr) - both of which are not 
very illuminating and we do not explicitly write them out here. Finally, we choose the 
boundary condition 

" d^~ 

(3.20) 



d u 



dT 



at the asymptotic boundary to account for the temperature derivative of the dimensionless 
mass parameter - giving the problem numerically slightly non-trivial mixed boundary 
conditions. Alternatively, the heat capacity can also be written in a more systematic 
way, 



B d p d 

c v = -T[d T + 2-— + 2^ — 

T dp 



AB 



2XB 



TdB 

„2 



M q d 

~f~IM n 



eT 



m M ( d d „ r d 

8pB— + — 4B— + 4p- M Q -r— 

T T \ dB ' dp q dM q 



C - 6 



(3.21) 



T ' 



where we identified the magnetic susceptibility, density of states, and mean magnetic mo- 
ment: 



XB 



dM \ 

dB \T,V,p 



-1 _ djj, 
£ — dp 



(vrT)4 / 

i67r — r U- 

■2 T 4 J \ °r 



dlC e + d § E*™ 



t,v,b ' \N 2 N 2 f T A 



m 



dM 
dp 



t,v,b T/\N c NfT 4 



dgd- p C e + dpE^, 



d^ 
dB 



(3.22) 
(3.23) 
(3.24) 



The derivatives of the scalar ^£ and ^| can be computed in an equivalent fashion as the 
derivative ^ described above. 
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3.3 Massless case 



In the massless case, the free energy can be integrated straightforwardly analytically and 
becomes 



F = — 



47rT 5 rg / du 
o 



u A -3 



+ 



- fp 2 + B 2 )u^ 



u 1 



V 



-V\T 3 N c N f ^ ( V /J + p + B~ 



(3.25) 



+ 2 (- (f + p 2 + B 2 )f" T (sinh- 1 (- (j 2 + p 2 + B 2 )) 



1/4 



- 1 



where the first term in the integral cancels the divergence of the second term at u 0; we 
substituted field theory quantities in the dimensionful factor in the second line and J-(-\-) 
is the incomplete elliptic integral of the first kind. For convenience, we give the asymptotic 
expansion F (sinh" 1 (-X) 1/4 \ - lj = i/C(2) - (-X)" 1 / 4 + 0(X)- 5 / 4 and the expansion 
at 



small values T (sinh- 1 (-X) 1/4 | - l) = 

(-X) 1 / 4 (1 - ro) + 0{-X) g / 4 , where /C(.) 
is the complete elliptic integral of the first 
kind and /C(2) ~ 1.854e i7r / 4 . These asymp- 
totic approximations are indicated in fig. 2. 
We also note again the effect of the electric- 
magnetic duality, which relates quantities 
under the interchange of density and mag- 
netic field. Now this emerges in the form 
that the free energy depends only on the 
variable f 2 +p 2 +B 2 . The representation of 
this symmetry in the various response func- 
tions then follows straightforwardly. For 
example the chemical potential is then re- 
lated to the magnetization as p(p,B,f) = 
only either of them. 



-0.5 



2 

— - -(-l)" 4 /f(2)/ 3/2 -f 
3 




1 f 




3 2 




F 









Figure 2: The free energy density of the defect 

as a function of \J f 2 + B 2 + p 2 (denoted for sim- 
plicity as "/". 

M(B, p, /), and in the following we will show 
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Now, we can simply verify the relations (3.10) and (3.11): 

1 



M = P 



J-^sinh- 1 (- (/ 2 + p 2 +B 2 ^ V4 



(- (p+p + B 2 )) 



1/4 



dF 



dF 



dF 



T irNcNjy/XT* Of 



= V\T 3 N c N f [ duA' t (3.26) 
Jo 

(3.27) 



1 F \sin\T 1 (- (f 2 +p 2 + B< 



1/4 



1/4 



= ttT 



/ dux ■ 
Jo 



(- (^p + p + B 2 )) 

This is shown in fig. 3, and we note that the dominant dependence of a on p or M on 




Figure 3: The baryon number chemical po- 



Figure 4: The contribution of the defect to 
the entropy as a function of \J p 2 + B 2 and 



tential function of p and y B 2 + p. 

Note that, due to electromagnetic duality, 
this is the same as the magnetization as a 
function of B and \J p 2 + p, upon an appro- 
priate scaling with a dimcnsionful constant. 



B, with an apparent saturation behavior as the scaling power changes from 1 to 1/2. The 
dependence on the subleading quantities is a small suppression, which goes slightly against 
classical intuition. Further more, we can straightforwardly give the contribution of the 
defect to the entropy, which is 



S = V\F 2 N c N f 



( 



\ 



jYsinh- 1 (-(f+p 2 +B< 
l+f 2 +p 2 +B 2 + 2f 2 ^- 



1/4 



P + p 2 + B 2 )) 



1/4 



/ 



(3.28) 

and is shown in fig. 4. Now, we can notice that for large values of /, this expression turns 
negative. Since the thermodynamic quantities that we derive here are only contributions 
to the overall quantities of the system "defect+background" , and the defect is certainly 
strongly coupled to the background and hence always in thermal equilibrium, this is not 
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troublesome, as we can demonstrate: 

First, we assume that the defect effectively extends roughly homogeneously over a region 
up to z max in the normal direction. This is suggested by (2.10), (3.26) and at large / also 
by the quasiparticle spectrum found in [23, 24]. Then, we consider the entropy density 
over this region of (3 + l)-volume, 



Saz 



( 



2\f\-KT*N c N f 



f z + p 2 + B 



1/4 



1 + f 2 + p 2 + B 2 



f.P ( sinh 1 (- (/2 + p + 



1/4 



- 1 



2/ 



(3.29) 



where we look in particular at the density of the negative term, -2TrV\T 3 N c N f f. In our 
limit of large N c and N c ^$> f 2 , this term is precisely (minus) the contribution q 5s s s £ M ■ This 
can be interpreted simply in the way that the extra degrees of freedom due to changing 
N c —7- N c + 5N C become fully available only after the brane falls into the horizon, and the 
positive term in the entropy of the defect describes the extra degrees of freedom contributed 
by the defect inside that region. Furthermore, there is no region in which the total entropy 
of the combined system is negative. 

At large densities or magnetic fields, or low temperatures, the contribution to the 



entropy is just 4irp or 



'XN c NfB 



respectively - indicating that the number of degrees of 



freedom is independent of the temperature, and proportional to the number of quarks or 
magnetic states. 

We can also straightforwardly compute the heat capacity: 



c v = 2\f\T 2 N c N f 



pT ^sinh 



~ \ \ 1/4 



1 + 2/ 2 



P + p 2 + B 2 



P + p 2 + B 



(- (^p + p + B 2 )) 



1/4 



P + p 2 + B 2 



(3.30) 



which we show in figure 5. Again, we notice that the heat capacity is negative for large /, 
but as in the case with the entropy, this does not signify an instability, as the contribution 
is much smaller than the heat capacity of the background SYM in the appropriate region 
around the defect. Also here, the negative contribution to the heat capacity simply indicates 
that the additional degrees of freedom only turn on gradually over some region ~ z max away 
from the defect. At small /, the heat capacity is just cy = 2\f\TN c Nf-j==^ + 0{p). 

This is somewhat counter- intuitive, as one normally expects an increasing heat capacity 
with increasing density but as the system is strongly coupled, we can only speculate about 
explanations for this behavior. We have to keep in mind that p is only the net density and 
there is always a finite density of quarks and gluons. 
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Figure 5: The contribution of the defect to 

the specific heat as a function of \J ' p 2 + B 2 
and /. 



Figure 6: The magnetic susceptibility of the 
defect as a function of B and \/ p 2 + f 2 . This 
is the same as the density of states as a func- 
tion of p and J B 2 + f 2 . 



Finally, we can, for example, look at magnetic the susceptibility, 



XB 



V\N C N 



f 



B 2 



2vr 2 T 



(3.31) 



+ 



B 2 



F ( sinh- 1 (- (f 2 + p 2 + B< 



1/4 



p 2 + B 2 + f 2 



f 2 + p 2 + B 2 



1/4 



This expression is always less than Xe? = \^\N c NfK 2 T 3 - the susceptibility of the defect 



without any of the parameters turned on as shown in fig. 6 - and behaves asymptoti- 

X]3 ^~^J^^ or Xb ^~^ /^^ ^ respectively. Hence, the defect is 



2VB 



cally as x W ( 1 L^gl , 



diamagnetic. Again, we can blame this on the strong coupling, which may increase the 
spin-spin interactions. Similarly, strong coupling may have the effect of suppressing the 
energy density of the plasma as we increase the quark density. 

3.4 Massive case 
3.4.1 Phases 

Now, let us look at the phase diagram in fig. 7. Looking at the different kinds of embed- 
dings, we find that in the case of vanishing p and /, there are three phases: the stable 

( BH) 

phases of the blackhole embedding at small masses M q < Mfnax (or large temperatures) , 
denoted as "B" and the Minkowski embedding at large masses M q > M^f 1 (M). In prin- 
ciple, there exists also an unstable phase with Mmax 1 > M q > (C), but in practice, 
this phase is not realized and there will just be a first order phase transition from B to 
M. At finite / or p, the mass diverges at u max = or = 1> so we find two additional 
phases, each one for the blackhole (Bl) and Minkowski (Ml) embeddings. We interpret 
this Bl embedding as a continuous deformation of the Minkowski phase with free quarks, 
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but also a high density of mesons. In many situations, the phase transition between those 
phases will disappear, as we will see below. Using the condensate as an order parameter to 
identify different phases, we show the various phases in figs. 7-9. In order to demonstrate 
the phase transitions, and the disappearance of the blackhole phase M and the transition 
phase C in the presence of strong magnetic fields, we plot the results at fixed temperature 
with varying M, rather than fixed mass and varying T. Obviously, in a thermodynamic 
process Mj (ttT) = cannot be attained. In fig. 7, we find that the critical mass decreases 
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The light gray (upper) surface is the mesonic, i.e. 
Minkowski, phase and the darker (lower) surface 
the blackhole phase. The gray line is the pro- 
jection of the minimum mass of the Minkowski 
phase M^l^\ The black line is the minimum 
mass of the blackhole phase and the dark gray 
line the local maximum Mmax 1 , if it exists, or 
the slowest rate of change of the mass with re- 
spect to the embedding. There is some numerical 
noise visible in the plots that can be ignored, but 
we chose not to suppress artificially. 

Figure 7: Value of the condensate as a function of the dimensionless mass, rh and magnetic field, 
B at fixed temperature. Top left: Vanishing density and vanishing compact flux /. The light and 
dark gray lines identify the transition between B and C and M, respectively. Top right: Density 
p = 0.5. The dark gray line is now the transition between B and C, or B and Bl and the black 
line between C and Bl. The light gray line is the transition between M and Ml. Bottom left: 
/ = 0.25. More details are explained in the bottom right. 



approximately linearly with the increasing magnetic field and there is a critical magnetic 
field, above which the blackhole phase B disappears. This indicates that the magnetic field 
catalyzes meson formation. This behavior is similar to what was observed in 3+1 dimen- 
sional systems in [32]. Rigorously speaking, however, there exists at all magnetic fields a 
continuation of the B phase at (exactly) zero mass and we will see in the next section that 
we can always attain both the Minkowski and blackhole phases at fixed magnetic field and 
mass and varying temperature. Also, we find how the phase C, that is suppressed and 
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even disappears at finite density or / re-appears at large magnetic fields. The first order 
phase transition from B to Bl turns into a smooth crossover as C disappears, and one 
might wonder what happens at the particular point, where C disappears. 

The condensate is large in the blackhole phase Bl and in the Minkowski phases and 
grows approximately proportional to the magnetic field - indicating chiral symmetry break- 
ing. As indicated above, these phases have a very similar behavior. The reader is reminded 
however that at finite density and vanishing /, the Minkowski embedding is not physical 
within string theory. It is interesting though that turning on / reverses the linear mass 
dependence of the condensate when interpreted as the conjugate potential of the mass - 
meaning that if we increase the mass or lower the temperature, the free energy will "sat- 
urate". In the blackhole phase B in which the mesons are dissociated and we have only 
free quarks, the condensate approximately vanishes and becomes independent of mass and 
magnetic field. 




Figure 8: Value of the condensate as a function of the dimensionless mass and density, at fixed 
temperature. Left: B = 0, right B = 4. 

In figures 8 and 9, we see first of all that the phase diagram at finite density and 
vanishing magnetic field is identical to the one at finite / under identifying p «->■ /. This can 
be seen explicitly also from the equation of motion (2.11) of the scalar, which is identical 
under the exchange p <H> / at B = 0. This is however only accidental and disappears 
as soon as we look e.g. at the value of the action or the free energy and the response 
functions. Turning on a magnetic field breaks this symmetry even though there is some 
overall similarity. It has the expected effect of lowering the critical masses, increasing phase 
space volume of the phases C and B and shifting the overall condensate - breaking also 
the symmetry of the M phase. Interestingly, this persists even at relatively large densities. 
To understand this behavior we recall that there is always a finite total baryon density, 
and the polarizing an localizing effect of the magnetic field is not affected by changing the 
ratio of quarks and anti-quarks. 
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Figure 9: Condensate as a function of the dimensionless mass and /, at fixed temperature. Left: 
B = 0, right B = 4. 

3.4.2 Stability and chemical potentials 

In order to determine which phases are thermodynamically preferred, we can look at the 
entropy. Before studying the data, we should notice that there is a trade-off between 
numerical accuracy and noisiness. This arises because a significant contribution to the 
entropy comes from the UV regime, i.e. from small values of u. On the other hand, the 
solution for \f becomes unstable and noisy in this regime, so we usually choose a cutoff 
Umin of the order of 10 -5 , which causes usually no significant errors in the result - except 
for the case of the entropy. Even if we try to extrapolate at small u, the cancellation of the 
boundary term in (3.15) will not be accurate, so we need to push the minimum value for u 
as far as possible and we will notice some noise. Eventually, the qualitative result will not 
be affected in either case, and we can further check whether some apparent "effect" is due 
to numerical errors or not by tuning Umin m 

We now use quantities that are made dimensionless using the mass as we consider 
the system at fixed mass. For example, we have T = i and B = -=^. We chose this 
combination, because these are parameters that naturally arise in the computations. Notice 
however, that m contains a factor of V% he. T = V~^ 2 3/2 M and B = = X g7T 2 M i ■ 

In the results for the entropy in fig. 10, let us first look at a few technical issues. In 
the plot on the top right, at / = 1/2, we notice that there is some numerical error which 
causes the entropy of the lower branch of the Minkowski embedding to remain finite at 
vanishing temperature. In the figure on the bottom-left, we "cut out" the region of small 
temperatures and densities, as it was dominated by noise because the configuration in this 
region requires embeddings with ~ 1 5 which are numerically problematic. Also, in the 
plot in the bottom right, the hint of the entropy of the phases crossing around T ~ 0.25 
can be shown to be due to numerical errors. 

Looking at the results, we see that in all cases, the blackhole embedding is preferred. 
In the case of having only the magnetic field non- vanishing, we actually notice the cusp-like 
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Figure 10: The dimcnsionlcss entropy S — 8M 2% Nf as a function of T = \ZA 23/ T m and B = 
^ %7r ? M 2 for / = = p (top left) and / = 1/2, p = (top right) or as a function of T and p (bottom 
left) or / (bottom right). 

behavior at the point where the blackhole and Minkowski phases meet that was already 
observed for (3+1) dimensional systems in [8]. We also notice again in the bottom right, 
which shows the plot "from behind", that / causes a negative entropy contribution, as 
discussed in the context of fig. 4. The fact that the difference in entropy between the 
Minkowski and black hole phases vanishes as T — > in the case of having the internal flux 
/ turned on, but that it doesn't vanish for the embeddings at finite density indicates that 
the latter are "more unstable" . This reflects precisely the observation that Minkowski em- 
beddings at finite density cannot be supported in string theory - in contrast to Minkowski 
embeddings at finite flux /. From another point of view, the vanishing entropy at T — > is 
what we expect generically and the finite entropy implies some "degenerate ground state" . 

Furthermore, we can see how the phase diagram changes if we consider a fixed quark 
mass, rather than a fixed temperature, m — > gets mapped to T — > oo and is hence 
not accessible anymore, and in fact finite values of B also get mapped to B — > oo in that 
limit. Hence, at all finite values of B, there are both blackhole and Minkowski embeddings 
and the critical temperature T c of the phase transition is some monotonously increasing 
function of B. Since there is no indication of a transition between the phases B and Bl 
in the plots in figure 10, let us see why. In our new variables a trajectory of fixed B and 
varying T corresponds to a trajectory B = A in fig. 7. Hence we need to look at e.g. 
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Figure 11: Left: "Zoom" to the region of the phase transition between the different blackhole 
phase embeddings Bl (low temperatures) and B (high temperatures) at small values of the density 
and fixed magnetic field B = 2. Right: the phase transition at large magnetic fields and fixed 
density p = 0.5 



large values of B and T at fixed p or / and the critical temperature will be T c oc Z? 1 / 2 . If 
we expand the equation of motion for \I> around small values of ^ and study the resulting 
linear second order equation, we find that this phase transition always exists for sufficiently 
large B and T. For finite values of / > 1, the onset of the phase transition gets shifted 
to B,T > 1, whereas for finite p > 1, the temperature scaling of p implies that the phase 
transition always appears at finite values of the temperature. In fig. 11, we show how the 
phase transition appears at small values of p and fixed B = 2, and at increasing values of 
B for fixed p = 0.5. 

In the light of the fact that the blackhole embedding is always thermodynamically 
preferred, the phase Bl and the phase transition between B and Bl are a smooth contin- 
uation of the Minkowski phase and the blackhole- Minkowski embedding phase transition 
as we we turn on the quark density p. This is also reflected in the fact that the Bl phase 
has lower entropy due to the larger fraction of quarks bound in mesons. 




Figure 12: The magnetization M as a function of the temperature-mass ratio T and the magnetic 
field B at p = = / (left) and p = 0, / = 0.5 (right). 
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Let us finally just take a quick look at the other derivatives of the free energy. In 
fig. 12, we show the magnetization and we notice that the magnetization is higher in the 
blackhole phase, which is dominated by free quarks, and lower in the Minkowski phase 
that are dominated mesons. In the case of small finite / (qualitatively the same happens 
also for p) however, the difference is highly suppressed and both phases have essentially 
the same magnetization. 




Figure 13: The chemical potential /i as a function of the temperature-mass ratio T and the density 
p at B = = / (left) and as a function of B at p = 1/2 (right). 

In figure 13, we look at the chemical potential. In contrast to the magnetization, there 
is a a significant difference in the chemical potential, that persists in all cases, with the 
one of the Minkowski phase being higher than the one of the black hole phase. Hence, 
inducing a difference in the density of quarks and anti-quarks requires more energy in the 
Minkowski phase. This arises intuitively, as the Minkowski phase is dominated by mesons. 
It might also be related to the fact that in string theory, there should be no finite net 
baryon density in the Minkowski phase. 

In the remainder of this section, we look at the second derivatives of the free energy 
that we obtained with the methods described in section 3.2. In figure 3.4.2, we show the 
magnetic susceptibility in the presence of the various parameters. In each of the plots, as 
also in figures 3.4.2 and 3.4.2, we show also the / = p = B = case for reference. In 
the top left, at finite magnetic field, we observe the saturation behavior with increasing 
magnetic field and the phase transition with the temperature independent behavior in the 
Minkowski phase and the "dimensional" T~ l scaling in the blackhole phase. This is because 
the Minkowski phase is dominated by the physics and energies of mesons, whereas the 
blackhole phase is dominated by the statistical mechanics of free quarks in the plasma. Also, 
turning on a small density in the presence of the magnetic field changes the susceptibility 
only sightly, in particular around the region of the phase transition, but turning a finite / 
changes it more significantly. Hence, apparently the net quark density is less relevant than 
the total quark density. Also the large temperature "tail" is universal, as it is dominated 
by the thermal equilibrium of "quark-anti-quark" production. At vanishing magnetic field, 
we see no divergent susceptibility at the phase transition, which is expected as there is no 
remnant magnetization at B = and hence also the "latent magnetization" of the phase 
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transition vanishes. In the presence of (only) a finite density, it appears that the high 
temperature blackhole phase also extends to small temperatures. For numerical reasons 
as we will discuss in the below, we cannot trace this phase to vanishing temperatures. 
However, the one of the fiducial Minkowski phases that typically mimics the blackhole 
phase is finite at vanishing temperature. As we saw in the massless case in section 3.3, 
increasing the density lowers the magnetic response. In the presence of (only) finite /, 
the behavior is significantly different. On the one hand, the curves do not converge in 
the high temperature regime, mainly because of the absence of a temperature scaling of 
the parameter /. At large temperatures, we still have the effect that finite / lowers the 
susceptibility, but at small temperatures, the susceptibility increases. Tracking the "lower" 
Minkowski phase to T — > shows actually that the lines of / = 1 and / = 4 cross and 
the zero temperature magnetization remains finite. The (stable and hence more relevant) 
blackhole phase however appears to diverge at small temperatures - but the numerics start 
to fail in this branch around T ~ 0.05 (4^0 ~l-5x 10 -4 ). This means that the defect is 
becoming magnetized. However, our methods are not sufficient to explore whether there 
is some spontaneous magnetization at vanishing temperature. 

The "density of states" , i.e the inverse of the response of the chemical potential to the 
net quark density, is shown in figure 3.4.2. In the Minkowski phase, e is highly suppressed 
at some small, but constant, value as we expect for the mesonic phase. At finite magnetic 
fields, this increases with some multiplicative factor. This just means that the chemical 
potential decreases with increasing magnetic field, and there is some rapid increase as 
we turn on a small density. At large temperatures, however, all curves converge to the 
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/ = p = B = case. As we turn on a finite density, this suppressed density of states is 
lifted to larger, but still approximately constant values, now in the blackhole phase that 
"allows" free quarks. That is the case if inducing a net quark density corresponds to 
increasing the total quark density - with a finite quark mass. If we turn on a magnetic 
field in addition to the finite density, there is a strong increase in e in the low temperature 
blackhole phase. A wild speculation would be that this is related to the density of Landau 
levels. Just as in the case of the susceptibility, the high-temperature behavior depends 
significantly on /. This may be simply because increasing / increases the physical "width" 
of the defect. As it is the case with the other quantities, there is no apparent phase 
transition in the blackhole phase. Now it seems, that in the blackhole phase e vanishes 
as T — > 0, which is similar to the diverging magnetic susceptibility. In the presence of 
a small magnetic field, surprisingly, the low temperature behavior does not get scaled by 
some factor as in the case of / = 0, but it actually gets substantially shifted to some larger 
value. 

Finally, we study the results for the heat capacity in figure 3.4.2. There, we see that 
there is no significant difference between the Minkowski and blackhole phases, except for 
the diverging heat capacity around the phase transition and some shift between the phases 
in the / = p = B = case with a higher value in the blackhole phase, because of the extra 
degrees of freedom of free quarks. Increasing the magnetic field slightly lowers the heat 
capacity, as does turning on a finite density. We can also nicely observe the peak in the heat 
capacity around the crossover between the high and low temperature blackhole phases at 
p = 1, B = 4. There are a divergence at T = and some noise at small temperatures, both 
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of which are due to numerical inaccuracy. There are two more noteworthy observations. 
One is the negative heat capacity at finite /. This effect is simply due to the fact that the 
extra degrees of freedom from modifying N c — > N c + 5N C are not available from z = but 
turn on only gradually over the defect - as we discussed in the context of the entropy in 
section 3.3. The numerics at / = 1, B = 4 at small temperatures seem to be not reliable. 

The other interesting observation is the fact that the heat capacity vanishes very 
quickly with decreasing temperatures in the blackhole phase at finite densities. As men- 
tioned before, the numerics let us down in this case at T ~ 0.7 in the p = 4 case and 
around T ~ 0.74 for p = 1, but already at these values c v has converged to much faster 
than expected for some power law scaling. This region is very interesting, as it turns out 
that at fixed p = J^, large values of rh = ^ are obtained at some "finite" G]0, 1[ and 
not at — > 1- This arises because at large m and p, we have rh(^o,p) oc p 2 . The value 
of determines the proportionality factor, such that we get a value of depending on 
p. Such coincidence in the scaling arises generically only in dd e f. + 1 dimensional matter 
embedded in dbuik + 1 = 2^/. dimensions. From a physical point of view in the gravity 
picture, the probe brane becomes "stiffer" with increasing p, which causes some "straight" 
limiting profile as we take T — > 0. This limit can however not be suitably studied with 
the methods of this paper, as we tune (or u max ) to obtain rh. A similar behavior also 
occurs in the corresponding branch of the (fiducial) Minkowski embedding. 



4. Discussion and Conclusions 

In this paper, we applied the holographic duality to study the thermodynamic properties 
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of a class (2 + l)-dimensional defect CFT's emersed in a heat bath of 3+1 dimensional 
M = 4 SU (N c ) SYM theory. To do that, we considered the properties that are extrinsic in 
the "area" of the defect, i.e. from a (3 + l)-dimensional point of view are independent of 
the placement of boundaries at some z^y away from the defect - in the limit of large z^y- 
We considered primarily the string theory setup realized by embedding N { probe D5- 
branes in an AdSs x S 5 background, in which case the system (at T = 0) preserves eight 
supersymmetries. At vanishing quark mass, the thermodynamic properties are identical 
(modulo overall factors) to the defect CFT dual to embedding N t probe D7-branes in the 
AdSs x 5 5 , which preserves no supersymmetries. As this contains only fermions at the 
massless level, it is of particular interest to condensed matter physics. As discussed in [24] 
the results of these embeddings are not reliable because it is not quite sure in how far the 
D7 setup persists in the light of gravitational backreaction. 

As we outlined in section 2, we considered deformations to the theory corresponding to 
a shift N c — > N c + 5 N c in the level of the gauge group on one side of the defect parametrized 
by/ 

~~ VJNf ' to a fi m * e quark mass M q , finite net baryon number density p and magnetic 
field B. 

First, in section 3.1, we discussed the choice of thermodynamic variables. We Legendre- 
transformed the action in order to use the density as a thermodynamic variable, rather than 
the chemical potential. There we also demonstrated that the parameter / is represented by 
the change in (3+1) free energy and its dual, defect width. As / also parametrizes the shift 
in the level of the gauge group SN C we chose to keep it as a fixed parameter, for which we 
also needed to Legendre transform with respect to d u z. In that context, we demonstrated 
how the tension of the probe brane at its endpoint in the case of Minkowski embeddings 
is matched by the appropriate number of attached D3 branes. 

We then derived the response functions, which allowed us to verify the expressions for 
the chemical potential and dual condensate assumed in section 2.2. We also developed a 
method to obtain the second derivative of the thermodynamic potential in closed-form, i.e. 
without taking numerical derivatives. 

In the massless case, we were able to obtain analytic expressions for all the thermody- 
namic quantities. We presented those in section 3.3, were we also discussed the negative 
entropy of the defect, and verified that it does not imply a negative entropy region and 
is simply related to the fact that the extra degrees of freedom due to 6N C turn on only 
gradually over the "width" of the defect. We also observed the implications of the electro- 
magnetic duality that was discussed in [24]. Here it corresponds again to an exchange of 
dimensionless density and magnetic field, p o B, and e.g. an exchange of chemical poten- 
tial and magnetization. Overall, we found a saturating behavior, where e.g. the chemical 
potential scales as p oc yfp at large densities. 

In section 3.4, we discussed the massive case. Using the dual condensate as an order 
parameter, we mapped out the phase diagram. We found that c vanishes in the blackhole 
phase but not in the Minkowski phase, where it increases with increasing magnetic fields; 
and that the critical mass of the phase transition decreases, i.e. the critical temperature 
increases, with increasing magnetic field similar to what was observed in [32, 33]. In the 
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case of a finite density or finite /, the blackhole phase extends to large dimensionless 
masses, i.e. vanishing temperature and the Minkowski phase splits into two metastable 
phases that extend to large mass or vanishing temperature. This deformation is smooth, 
and at small values of p and / the blackhole phase tracks closely the Minkowski phase, 
and hence we find a transition between two different blackhole phases. At larger p, /, 
this transition disappears and also the two Minkowski phases split further - one remains 
approximately unchanged and the other one approaches the blackhole phase. The phase 
transition, or minimum mass/largest temperature of the Minkowski phase shifts towards 
lager masses, approximately fh C rit. oc \fp, \/7- 

We then studied the results for the response functions. To do so, we considered fixed 
mass, and quantities made dimensionless with the appropriate powers of the combination 
2 a/2/ XM q , rather than ttT. First, we looked at the entropy and found that the blackhole 
phase has higher entropy than the Minkowski phase, as expected since it is dominated by 
free quarks rather than mesons in the Minkowski phase. The entropy vanishes at vanishing 
temperature within the limitations of numerical accuracy, except for the blackhole phase 
at finite density, where there is some "ground state degeneracy". We also demonstrated 
how the phase diagram gets mapped in this parametrization of fixed mass; e.g. the phase 
transition at large magnetic fields is at T oc y/~B. 

Overall, in the high temperature regime, the effect of fixed B and p disappears, because 
the relevant parameters are the dimensionless B and p that scale oc T~ 2 , but the effect of 
fixed / obviously remains. The behavior in this regime is verified to be consistent with the 
massless limit. 

At small temperatures, there are a few interesting effects. At finite / the crossover 
between the blackhole phases disappears and one of the Minkowski phases is very similar to 
the blackhole phase. For the magnetic response, we find that at finite /, the defect seems 
to become magnetized at vanishing temperature, as the susceptibility diverges at small 
temperatures and vanishing magnetic field but is finite and constant at finite magnetic 
fields. In the Minkowski phase, it becomes large, but finite. A similar behavior occurs 
with the response to the density, where the "density of states" , the inverse of the second 
derivative of the free energy, vanishes. Certainly, it would be interesting to see whether 
there is a spontaneous symmetry breaking at vanishing temperature and finite M q and /, 
but this is beyond the scope of this paper. 

Atp = B = f = and in the presence of finite magnetic fields, the density of states 
below the phase transition is small but finite and approximately temperature independent 
and increases with increasing magnetic field. The susceptibility is also approximately 
constant, but of 0(1) and decreases consistent with the magnetic saturation effects. In the 
case of finite density, the density of states is also 0(1) and approximately constant at low 
temperatures, increasing with increasing density consistent with saturation effects and the 
susceptibility is finite and increases with decreasing temperature in a smooth continuation 
of the high temperature behavior. 

A big surprise however happens for the heat capacity, which vanishes at finite p with 
decreasing temperatures more quickly and at larger temperatures than expected for the 
general T 2 behavior. This is particularly interesting in the context of the unusual embed- 
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ding in this regime. Rather than at — > T i- e - the limit that the probe brane intersects 
the horizon only at a point, we obtain T — > at fixed M q and p at some smaller value of 
*o- This arises because of a coincidence in the dimensionless scaling of the quark mass 
M q /T oc \J p/T 2 at small temperatures, with a proportionality constant depending on 
and hence relating p/M 2 to ^o- Physically, this means that the brane has a higher effective 
tension, i.e. becomes "stiffer" as we decrease the temperature at fixed p, i.e. as we increase 
p/T 2 in such a way that there is a limiting "straight" embedding at fixed as we decrease 
the temperature at fixed mass, i.e. as we increase M q /T. It would be very interesting to 
explore this new behavior further, but we will postpone this to future research as the meth- 
ods used for this paper are not suitable and break down at small, but finite temperatures 
in this regime. 
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